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ON THE SELECTION OF THE CLASSICAL LIMIT FOR 
POTENTIALS WITH BV DERIVATIVES 



AGISSILAOS ATHANASSOULIS AND THIERRY PAUL 



Abstract. We consider the classical limit of the quantum evolution, with 
some rough potential, of wave packets concentrated near singular trajectories 
of the underlying dynamics. We prove that under appropriate conditions, even 
in the case of BV vector fields, the correct classical limit can be selected. 
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1. Introduction 

The fundamental equation of quantum mechanics is formulated either as the 
Schrodinger equation for a wavefunction, 

(L1) ie-§ i uZ(t) = (-^A + V(x))u%t), 

u e {t = 0) = u§0), 

or more generally as the Heisenberg-von Neumann equation for a density matrix 



(1.2) 



iejj i D s (t)= -^A + V,D s (t) 
D%t = Q)=Dl. 



The state of the quantum system is described by the operator D £ (understood to be 
the orthonormal projector |u e )(it e | when working with the Schrodinger equation), 
evolving in time under the aforementioned equations. Indeed, under very general 
conditions on the potential V , known as Kato's conditions, this is a problem well 
posed for Uq G L 2 or Dq being a Hilbert-Schmidt operator HH US] ■ It is also well 
known that the positivity and trace of Dq is preserved in time. The parameter e is 
called Planck's constant, and when it is very small one usually expects the system 
to behave like a classical one. 



2 



ATHANASSOULIS AND PAUL 



An equivalent way to write equation (|1.2[) is in terms of the Wigner transform, 



(1.3) W £ (x,k,t) = W £ [D £ {t)]{x,k) 



y 



where for each time t K £ (x,y,t) is the integral kernel of the operator D £ (t); in 
other words, for / e S 

,x + y 



D £ f(x) = e 



: W £ {—^,k)f{y)dydk. 



A compact way to say that is that the Wigner transform of an operator is e" times 
its Weyl symbol. Thus the operator corresponding to the Wigner function Wq is 



l O PWeyl (W £ ). 



In the case of a pure state this yields 
(1.4) W £ {x,k,t) = W £ [u £ {t)]{x,k) = 



e - 2myk u s (x + e^,t)u £ (x - e|,i)dy. 



Equation (1.2 1 implies 10J 
(1-5) 

d t W £ + 2irk ■ d x W £ + i J e~ 2mS y v(x+ i y) : V ( x +W dy W £ (x, k-S, t)dS = 0, 
W £ {t = 0) = Wfi = W £ [Dq]. 



The propagator for equation ( 1.5 1 is constructed from the one of equation ( |1.2[ ) and 
as long as Kato's conditions apply and Dq is a Hilbert-Schmi dt op erator (equiv- 
alcntly Wq E L 2 ) there exists a unique solution for problem (1.5), and it is the 
Wigner transform of the solution of equation ( 1.2 1; see theorem 5.4 below, and [IT]. 

We will often use the shorthand Tj for the operator involving the potential (see 
Definitions and Notations in section [3] below) . 

The Wigner transform (WT) respects the structure of the problem, allowing the 
Wigner equation (1.5) to inherit properties (and in particular conservation laws) 
from (1.2), see e.g. lemma [5^4] Moreover, under appropriate conditions, it allows 
for a very natural and compact description of the semiclassical limit, e -> 0. Indeed 
the WT has a physically meaningful limit as e tends to zero, while in general the 
wavefunction u £ itself (or the operator D £ ) does not. The limit (in the weak-* 
topology of an appropriate algebra of test functions), called the Wigner measure, 

W £ W°, 

satisfies the Liouville equation of classical mechanics 

d t W° + 2nkd x W° - ±d x V ■ d k W° = 0, 

(1-6) 

W°(t = 0) = UmWfi 

if the potential is regular enough jTDl [7] . 

For potentials with low regularity (less than C ' ) two different problems natu- 
rally arise; 

Problem 1: showing that the Wigner measure is a weak solution of (1.6); 
and, if possible, 

Problem 2: constructing a selection principle that identifies the correct weak 
solution, since problem (1.6) can be ill-posed in that case. 
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Problem 1 is solved in [10| for V € C 1 (here we will work for quite less regular, 
if somewhat specific, potentials). As is highlighted already for V £ C 1 \ C 1,1 , there 
can be several weak solutions when the initial datum is a singular measure (i.e. 
supported on a set of measure zero, the most natural case being a single point). 
That is what we mean in the sequel when we say the quantum initial data Wq is 
concentrating - concentrating to a point-supported measure W§ = lim Wq. Indeed, 

£->0 

the Liouville equation can be made well-posed for much worse potentials than C , 
but typically not for concentrating initial data [3 [51 [HI E] • 

More recently, Problem 1 has been successfully handled in [S] under BV con- 
dition of the vector field generated by the classical Hamiltonian, and for non- 
concentrating (in the sense defined above) and small (in operator norm) initial 



data. In that situation the solution of (1.51 tends weakly to the push forward of the 



limit of the initial datum by the so-called DiPerna-Lions-Bouchut-Ambrosio flow 
(defined a.e.). In 0] a new technique was formulated for handling semiclassical 
limits with concentrating initial data. Parts of it were used to work out problems 
with low-regularity potentials (also in strong topology), in [3j; see also [12; for a 
short review. In a nutshell, the technique of [3] applies to more general initial 
data than the results of [31 [3] (including in particular data concentrating to a 
point-supported measure), but less general potentials (roughly speaking, of C 1 '" 
type). 

In the present note we compute explicitly the Wigner measures for certain prob- 
lems with potentials that are not in C 1 , but have measure valued second derivative 
(therefore in a sense having the typical flavour of BV vector fields), and initial data 
which concentrate on the points of singularity. We illustrate the idea behind this 
computation and discuss how this idea can be generalized to other problems with 
similar regularity, and, roughly speaking, isolated repulsive singularities. The main 
difference from [3] is that here we manage to work with other topologies, more 
appropriate for the problem, and this is crucial to achieving stronger results. 

The main result of this paper is theorem |2.1| let us also mention the simple 
observation of Lemma 15.31 



2. Main result 

Since the regularity condition on the potential does not insure unicity of the 



flow, we will compare the solution of the quantum problem (1.51 to the behaviour 



of the solution of the Liouville equation with the same, e-dependent, initial 



condition. We will show that the limit of the solution of the problem (1.5| is the 
same as the limit of the solution of the Liouville equation with the quantum (i.e. 
e-dependent) initial datum. In other words, the extra information needed for the 
selection of the weak solution that correctly captures the semiclassical limit is fully 
contained in the way in which the initial datum concentrates to a point-supported 
measure. 

Theorem 2.1. Let ip : [0, oo) -> [0, 1] be a C°° monotone cutoff function, i/j(s) = 1 
for all s < k) tp(s) = for all s > 1. Now, for x € R 2 , some 9 € [0, 1) set 

v(x) = (i-\x 1 \ 1 + 8 ) v(NMN). 

Moreover let 

-(0,-A) k-(0,KY 



s x 
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forsomew £ H 2 n L°° nW 1 ' 1 , w(x,k) > 0, / wdxdk = 1, supp w C {x 2 + k 2 ^ 1}, 
X > 1,K > and — S 2 — (—log^e*)) 7 ^ (see claim 4.1 for other possible 
scalings). 

Let W £ be the solution of 

d t W e (x, k, t) + 2nk ■ d x W e (x, k, t)+ 

+l J e -2«iSy V(x + Sv)-V(x+tv) dy fc _ ^ ^ = ^ 



(2.1) 



W e (t = 0) = F £ 



(Recall that W 6 is the WT of the solution of equation (1.2) with initial data 
e~ n Opweyi (Wq)). Denote also by p e the solution of 

d t p e + 2irk ■ d xf f ~ ^d x V(x) ■ d k p e = 0, 

(2.2) 

f(t = 0)=F§. 



(With this regularity of V problem (2.2) has a unique solution in LP [T]^. 

= £ AnL 2 . 



Then, for each t £ [0,T] 



(2.3) 



lim{W e {t)-p s {t), 



(See Definition 3.1 below for the algebra A). 



Ln other words the semiclassical limit of (2.1) behaves exactly in the same way 



as the concentration limit for the Liouville equation (2.2) 



It can be noted that the only kind of initial data allowed in theorem |2.1| are 
wavepackets concentrating on a given point. One reason the result is phrased the 



way it is, is that approximation (2.3) can hold even when there is no unique semi- 



classical limit. If, for example we restrict the values of £ £ { — } and we substitute 

— — j_ j_ _i 

the initial datum Fq — Ftf" by F m {x\ — d"> , X2, fci, fca) for dm = C(— l) m (logm) 2 , 

then it is easy to check that the result applies, but W™ (t) - and p™(t) with it - 

has two accumulation points, one scattered to the right of {x\ = 0} and the other 

to the left (see also next remark). The finding here is that, whatever the interaction 



with the singularity, it is the same for the quantum problem (2.1 ) and the classical 



problem (2.2). In other words, all the information needed to determine the interac- 



tion is contained in the initial datum. In that light, even this limited pool of initial 
data contains enough different possibilities to explore. 

We can use the theorem to compute in more detail the semiclassical limit: 
Corollary 2.2. 

(2.4) W°(t) = c + S(X+(t),V + (t)) + c_S(X-(t),T-(t)), 

with c± — J w(x, k)dxdk and (X ± (t),V ± (t)) are the two diverging curves ob- 

±x 1 >0 

tained by (X ± (t),V ± {t)) = lim(X±(t),V±(t)), 

X±(t) = 2*7>±(t), £±(t) = -&W(X±(t)), 
tf{t = 0) = (±n,-X), V±(t = 0) = (0,K). 
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This can in fact be used to construct quite exotic examples; in figure [T] one sees 
various trajectories (their projection on the X plane, to be precise) for various 
values of the parameter K. (The value of the parameter X is not really interesting, 
since it merely determines the length of free motion before interaction with the 
potential starts). Observe moreover, that by an appropriate superposition (discrete 
or continuous) of initial data as in theorem 2.1 one can construct an example where 
the outgoing waves are (discretely or continuously) distributed among almost any 
angle. 




Figure 1. Pairs of trajectories (X ± ,'P ± ) corresponding to var- 
ious values of the parameter K. The projection to the (X\^X-i) 
plane, as well as the contours of the potential V(X\,X<i) are shown. 



This is a highly exotic example, and it may well turn out that its physical rel- 
evance is limited. However, it does point out a couple of points that could prove 
rather fruitful. First of all, the tools of semiclassical analysis can be extended and 
applied even to phenomena that are qualitatively very different from the smooth 
regime that they originate in. Such splitting of a particle and other possible ex- 
otic examples can motivate technical (and possibly numerical) experiments and 
refinements on existing methods and tools. 

Moreover, this example points quite naturally to a more general (and, one ex- 
pects, physical) situation: there is a "scattering process" happening "on the sin- 
gularity" . That is, scattering with respect to local fast variables a neighbourhood 
of the singularity. In this particular case (of very slow concentration, since S x ,5k 
are logarithmic in the semiclassical parameter) the scattering operator is "trivially" 
given by the concentration limit of the respective classical problem. That is, it is an 
operator fully determined and constructed by classical dynamics. In cases of faster 
concentration (e.g. pure states) one expects that an auxiliary quantum problem will 
need to be solved in a neighbourhood of the singularity to determine the "splitting" . 
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In other words the irregular potential makes the corresponding flow to behave in 
O(l) time-scales as a regular one would in long times (diverging trajectories), and 
this comparison can yield some very insightful finds [13] . 

It should also be noted that this particular example rests on the exact alignment 
of the direction of propagation to the line of singularity (even a small rotation can, 
potentially, destroy it). As can be easily chec ked, if the potential is substituted 
by V(x) = (l — |a;| 1+6 ') Vd^iDV'da^l) theorem 2.1 still holds, and then is robust 
to any perturbation of the initial data. Still, the only nontrivial case is when 
the particle passes exactly over the singularity, and the possibility of non-standard 
effects in that case. 

At this point we must refer once again to [2]; a stable result, which covers "almost 
all initial data" - in which case the a.e. unique solution of the Liouville equation 
describes the semiclassical limit, and there is no room for exotic effects. In other 
words these interactions may be "atypical" in some sense, but if one tries to study 
them, they have to select initial data so that they fully interact. 

Let us finally remark that an outline of how to check whether this technique 
applies to different problems is given in section [6] 



3. Definitions and Notations 
The Fourier transform is defined as 

(3.1) f(k)=^ k [f(x)}= [ e- 2 * ikx f(x)dx. 



For compactness, we will use the following notations: 



(3.2) 



F 2 W{x, K) = j t 



-2-KikK 



W(x,k)dk. 



(3.3) 



TYW = ife 



-2-n-iSy v(,x+iy)-v(x+jy) 



dy W(x, k - S)dS 



= T 



K^k 



V {x+iK)—V (x— %K) 



T 2 W(x,K) 



= T 



K^k 



J VV(x + s)ds T 2 W{x,K) 



(3.4) TfW = ~^d x V ■ d k W = T K \ k [VV{x) ■ K F 2 W{x,K)] . 

Definition 3.1. The space A is defined as the completion of the smooth functions 
of compact support C^°(IR 2n ) under the norm 

II/IU = / sup\T 2 f(x,K)\dK. 

J X 

K 

It follows that 

II/IU' = sup f \F 2 f{x,K)\dx. 

K J 
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Some basic properties that we will use are 



< \\f\\A 

< WfWxA < \\f\\M 



Definition 3.2. Every finite signed measure can be decomposed to positive and 
negative part, = /i + — /i - , for some finite non-negative measures /x + , (JT . We will 
denote the total variation of a signed measure 



WiA\m = J + J dfi . 

MP R n 

Definition 3.3. We will denote by $ = $ £ the smoothing operator 

& : f(x) ^[IY J e-^ lJ ^f{x')dx'. 
When there is no danger of confusion, we will write 

f:=$ e f. 

Definition 3.4. The Sobolev space W m ' p is defined as the completion of the smooth 
functions of compact support C£°(R™) under the norm 

wm , P = iiv a /iu,. 

\A\^m 



4. Proof of the main result 

4.1. Strategy of the proof. The intuition behind this proof is quite simple: the 
quantum initial datum, while concentrating to a point still is an I? function, hence 
the a.e. theory for its evolution under the Liouville equation applies. The concen- 



tration limit of problem (2.2) is therefore a natural candidate for the semiclassical 
limit. Of course to show a semiclassical approximation, a certain degree of smooth- 
ness in the potential and the Wigner function is needed, and in any frontal approach 
to this problem, such smoothness simply is not there. 

A simple idea to try out is the following: can we cut-off the pieces of the Wigner 
function that approach too closely to the singularity? If we do that, can we meaning- 
fully use an auxiliary function supported just far enough away from the singularity 
so that it preserves enough smoothness itself - as well as allowing one to cut-off the 
potential's singularities? This was pretty much the program we followed before, in 
®. 

The new element here is that to strengthen that technique to potentials as bad 
as the ones we treat here, basically a much bigger piece of the initial datum would 
have to be cut-off, and we need to find a meaningful way in which such a cut-off 
introduces "small" errors. To do that, the positivity of the density matrix comes 
into play, and some very different considerations are needed. 

Claim [44] holds all the compromises that need to be made, and is the conclusion 
of a lot trial and error. If one assumes it and move on in a first reading, the flow of 
the rest proof should provide a reasonable motivation for why these computations 
have to be just so. 
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4.2. Proof of theorem |2.1[ To facilitate the presentation, let us introduce at 
this point a number of auxiliary functions (using the notations introduced in the 
previous section): 

(4.1) d t W e + 2nk ■ d x W e + T^W e = 0, 

W s (t = 0) = F e , 

(4.2) d t Wt + 27rfc • d x W{ + T^W{ = 0, 

W|(t = 0) = F{ = *F%, 

(4.3) d t W£ + 27rfc • d x W% + T^W£ = 0, 
W|(t = 0) = F| = $(l-^))Fo £ , 



(4.4) SjWI + 27rfc • a x W| + T/W| = 0, 
W!(t = 0) = Ft = (1 - iK^))F§, 

(4.5) d t( of + 27rfe • d xP \ + T^pI = 0, 

pl(t = 0)=FI, 



(4.6) + 2irk ■ d xP e + T^p e = 0, 

ff(t = 0) = F§. 



The function tfj is the same as in the statement of theorem 2.1, and R' will be set 
below. 

Obviously, we are going to use these functions as stepping stones, passing from 
one to the other with the appropriate topology each time. Because this topology 
cannot be always the same, the end result is formulated as is, in weak sense: for all 

(f> e l 2 n A 

| (w E (t) -p e (t),(t)) \ sc \{w s (t) - W! (t),4>)\ + \ (Wf (t) - wi (t) 
+\m(t) - Wi{t),4>)\ + \(w§(f) - Pl(t),4>)\ + \(p\{t)-f{t),4>)\ 

The point is to collect the various constraints that would come from each of these 
building-block problems and satisfy them at the same time. That is essentially 
done in the following 

Claim 4.1. With appropriate calibration of 5 x ,5k, R ', we have 

11^1-^11^=0(1), 

||pflk»([o,T],H 2 ) = o{e~i), 

and 

||*!-ff|| £a = (i). 

Proof of the claim: The point is to make sure that p\{t) doesn't pass through 
a neighbourhood of the set where the second and third derivatives of the potential 
are singular (or too large in any case); in this case the strip {|x| < R}, R = 
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(-log(e*))v=s , for t € [0,T]. In that case, using lemma |5.l| it follows readily that 



indeed 




(by recalling that in fact only the values of the potential's derivatives along the 
part of phase space that the solutions passes through matter). 

To ensure that p\{t) stays away from {|x| < R}, it turns out to be sufficient to 
exclude a somewhat larger strip {|x| < R'}. To see why, imagine firstly that we are 
in free space, V = 0. The largest possible momentum in the x\ direction is 5k\ in 
time T this can only cover a distance of TSk', therefore R' = R + TSk should do it. 



Figure 2. Any trajectory leaving the support of F§ has an initial 
velocity in a small cone around k = (0, K). If we were in free space, 
there wouldn't be enough time for sufficient movement in x\ to 
reach {|xi| < R} for t € [0, T]. Therefore the solution is the sum 
of two components, supported on either side of {|xi| < R}. The 
time-scale T enters in the constants used, but the construction can 
be made for any fixed T. 

Now we have to take into account the addition of the potential. The repulsive 
nature of the singularity (preserved by the cutoff ^(|xi|)) makes sure that any 
trajectory moving towards it will not be accelarted - but in fact slowed down; 
i.e. any movement towards the singularity will be less than TSk- The presence of 
V'(|x2|) serves only to slow somewhat movement in the X2 direction (or even turn 
it back); but makes no difference whatsoever in the X\ direction. 

The timescale T is interesting when trajectories starting in the support of Fq 
reach and leave the support of the potential; in any case T — 0(1). 

Moreover, for the first part of the claim, 





U \ 5, + S k + S x ' 
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Finally, for the third part, 

m-mv = w-* e )m*^vz\\m# = o^fe i +sj; i )%% 1 ) 



So collecting all the constraints, we have 

y/£<£6 x ,5 k ; R'<^S X 11^-^11^1=0(1), 

R' = (-log(si))<^ + C8 k ||pf |U~ ([0 ,t],^) = 0(e-i\\F§\\ B '), 

(S^ + S^S-^ 1 « e-i =* \\FSWn* = o(e~i), 

yfe « R'S 2 J k , y/i « R'S x 6l => - = o(l). 

A concrete scaling that makes all these constraints covalid is (recall that R 
(-Jos(ei))^) 

4 

The proof of the claim is complete. 



Sk = Si = R. 



Remark: The claim apparently depends on the geometry of the problem; but 
really all we used was that the flow is repulsive away from the singularity. That 
is, that any trajectory would not go towards the line {xi — 0} faster than it 
would on free space. Thus it makes no difference if e.g. a different potential - 
for which the same property holds - is used. Indeed if the potential was V(x) — 
(l — J ^(liEiD^da^l) nothing needs to be changed. 

Given this calibration, the proof proceeds in a very predictable fashion. Indeed, 



using the unitary propagation of the Wigner equation - theorem 5.4 - one observes 
that 

11*0-*? IU" < Ve\\F^\\ m = 0(VeS x S k (5 x +S k )) = o(l) \\W§{t)-W;(t)\\v = o(l) VI 

The estimate between W{ and Wf is one of the essential parts, and the techni- 
cal innovation here. Denote by U £ (t) the propagator of the Wigner equation; we 
observe that 

Wf-W§ = U e (t)(F?-Ff) = 

= um (f) 2 /e- 2 - '"' )2 ^ fc,)2 ^) d *' dk '> 

the point being that Ff — F| is itself a density matrix - convolution of a positive 
measure with a coherent state; see lemma |5.2| Moreover, it is a small density 
matrix, and this is preserved in time; see lemma |5.5| Indeed it follows by the 



conservation of trace (and lemma 5.3 for density matrices) that 

\W(t) - wnt)\u> = m - fi\\a- < m - f||Ui = o(y] 



which was already calibrated to be o(l) in the proof of claim 4.1 (Indeed following 



that calibration we have y- = 0{y/5 k + S x )). Observe that Ff — F| is not o(l) in 



5 

L 2 sense, hence the introduction of the L 1 -like norm is necessary. 
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Moreover, by virtue of theorem |5.4| and claim \4~T\ 

\\Wi(t) - W!(t)\\ L2 = \\FI - F§\\ L * = o(l). 

The other non-trivial step is the quantum-classical dynamics comparison between 
W3 and p\ . Recall that U e (t) the propagator of the Wigner equation, and denote 

by E(t) the propagator of the Liouville equation. Set h e (t) — W§ (t) — p\(t); then 

t 

of course h £ (t) = J U £ {t)(T^ - T^)p\(t)(1t. 

T = 

Using the Fourier transform in the k variable we have 

|| / U{t){TY -TX)pl{r)dr\\ L ^ 

T=0 

sup ||( g(g±gjzgfrzg) _ dxV (x) ■ K)F 2 pl(t)\\ L 2 < 
te[o,T] 

sup |[v -iMD C^^f^^ - d * V W ■ K)F 2 pm\\L* + 
te[o,T] lu,£ 1 

+T sup \\ x 1 (\K\)( ^±^n^fl _ dxV{x) . K )T 2P m\W < 

t£[0,T] (£ - +co > 



<e\\x, 0e - h (\K\)^ sup \d^V{x + r^)\F W m\\^+0{sh)\\\K\^ 2 pm\\^. 

- re (-1.1) 

\A\=2 



For the first term we used a straightforward Taylor expansion, while for the 
complementary case we used the computation 

|( ^ + fK(,-f) _ Wijf)|< 

, «|K1 

<l| / (^-t)J2z j z m d XjXm V(x + tz)dt\^C\K\j:\ } d x . Xm V(x + tz)dt\^ 

_ e\K\ j,m j,m _ e\k\ 

2 2 

<c\k\ E ll^llx. 

.4=2 

The key to proceed is to make use of the well prepared initial datum F% : denote 
S e = {x\3 t e [0, T],k G R 2 : p{(x,k,t)>0} 

and 

B e = (J - x*\ < y/e}. 

Then 

fK 

x ]Q - h {\K\) sup \ d £v( x + T— )K |^y(x)| 

l °' e 1 re (-1,1) 2 xeB? 

\A\ = 2 \A\ = 2 
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and therefore 



e X , 0e -h(\K\) sup \d£V{x + T<f)\ || \K\*T 2P m\\L- < 

lu ' £ 1 re (-1.1) 

|A| =2 

sCe sup \d£V{x)\ sup \\p 6 2 (t)\\ H 2 

x 6 B E te[0,T] 
\A\ = 2 



finally yielding 



/ ^(t)(T7-T^)p§(T)dr|| La 



r=0 



= 0(Vi||p|(*)ll^(l + Vi ^ |9^)l)) = 0( £ i||^|| H2 (l + J R e - 1 )) = (l). 

1-41 = 2 

This is possible precisely because our approximate initial datum F§ stays away 
from a neighbourhood of the singular set {\x\\ = 0} when pushed forward in time 
by the a.e. Ambrosio-Lions-Di Perna flow E(t), as was checked in claim 4.1 (For 
a potential of the form |a;| the second derivatives are zero almost everywhere, and 
this creates the possibility of taking advantage of a very special structure. That's 
why we included the case > in the statement of the theorem, to show that in 
principle this technique work for a variety of localized repulsive singularities). 

5. Auxiliary results 

Lemma 5.1 (2 nd order derivatives equations for the Liouville equation). Consider 
the Cauchy problem for the Liouville equation with potential V(x) on R™, 

d t f + 2nk ■ 8 x f - ±d x V(x) ■ d k f = 0, 

(5.1) 

/(* = 0) = /o- 

There are constants C\ , Ci > depending only on n such that 

tC 2 SUp \\d°V(x)\\ L oa 

\\m\\H> < Cie '»>« 3 H/oIIhs. 
Lemma 5.2 (Density matrices). Let p, be a probability measure on M. 2n . Then 

/2\ f n U-sc') 2 + (k-k')' 2 

W e (x, *)=(-) / e ' dp(x', k') 

is the Wigner function of a density matrix, i.e. the corresponding operator D e is a 
positive trace-class operator with tr(D E ) = 1. 

Proof: Though the proof is obvious by using Toplitz quantization, let us give a 
direct proof. 

We know that tr(D e ) = J W £ dxdk = 1. Let us now look at positivity. To that 
end, observe that the integral kernel 

K e ( X> y)=(~Y J e^^e'^-^+^-^dpix^ko) 

x ,k 
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is the kernel of a positive operator. Indeed: 



fK £ (x,y)u(x)u(y)dxdy=(Z)* J |( e - 2 ™^+? (*-*») , u )\ 2 d^x ,k ) > 0. 

The proof is complete by observing that the Wigner function corresponding to 
the kernel p £ is 

" -^V K e [x+ S JL yX _ £ JL )d y = W *. 

Lemma 5.3. For any trace class operator D e , with corresponding Wigner junction 
W £ , 

\\W e \\ A > < \\D £ \\ tr = tr(\D s \). 

Moreover, if D £ ^ 

\\W e \\ A , = \\D% r = tr(D e ). 

Proof. Let D £ be a non-negative trace-class operator. Then, it admits a SVD 
expansion over orthonormal projectors 

D £ = y^^X m \u m )(u m \, 

m 

where of course ||A m ||;i = ||Z? £ ||t r = tr(D £ ), i.e. A m ^ 0, and (u m ,ui) — 8 m j,. It 
follows that 

W £ = *>mW £ [u m ] = A m / e- 2 ^u m {x + £ f)u m (x - £ -^)d y 

m m 

and, by straightforward substitution, 

1 1 W £ | \ A , = SUp J | Z ^mUrn (X + £ ~f )Um (x - £ ~f ) \ dx < 

K x m 

m X a: m 



Now observe that 

||W e IU' = supj\Y,X m u rn {x+ E f)u m {x- E £)\dx> 

K x m 

> fY,^rn\u m (x)\ 2 dx = 
x rn m 

Hence, D e ^ || W e |U' = ir(£> £ ) = ||£> E || tT .. 
To conclude, observe that for D e = — Df_ 

||^ £ |U < \\W e [D £ + ]\\ A , + \\W £ [D £ _}\\ A , = tr{D%) + tr(Dl) = ||^|| tr . 
The proof is complete 

The following theorem is well known, and follows from the correspondence be- 
tween the density matrix and the Wigner function; see e.g. Theorem 2.1 of [IT]: 

Theorem 5.4 (L 2 regularity of the Wigner equation). If the Schrddinger operator 
— §A + V(x) is essentially self-adjoint on L 2 (WL n ), then the corresponding Wigner 
equation preserves the L? norm, i.e. for Wq G L 2 (M. 2n ) there is a unique solution 
of 

d t W e + 2-nk ■ d x W £ + ^Re 



(5.2) 

W £ {t = 0) = w§ 



i J e 2niSx V(S)W E (x, k- £ f)dS 



= 0, 
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and \\W £ (t)\\ L 2 {R 2 n) = \\W§\\l2(k^) Vt e R. 

Lemma 5.5 (Trace conservation). Let Wq be a Wigner function corresponding 
to a trace- class operator. Consider a potential V such that the corresponding 
Schrodinger operator is essentially self-adjoint. If by W £ (t) we denote its evolution 
in time under the corresponding Wigner equation, then 

tr(Wo) = tr(W e (t)) Vt. 

Proof. One easily checks that Wfi = £ X m \u m )(u m \ => = Yl X m \u m (t))(u m (t)\, 

rn m 

where u m (t) is the evolution in time of u m G L? under the Schrodinger equation. 
The result follows. 

6. Extensions 

As was mentioned earlier, this approach also applies if the potential is substi- 
tuted by V(x) = (l — V'd^iDV'd 2 ^!)- Other straightforward generalizations 
come by embedding the problem in higher-dimensional space (i.e. including more 
transverse dimensions) . 

Checking whether some version of this approach applies to a problem would start 
by building the counterpart of claim |4.1[ and its second part in particular. More 
specifically: 

• Set N = {(x, k)\ sup \d£V(x)\ > C} (C is a parameter of the order —loge 

|A|=3 

to some power). This is the "neighbourhood of the singularity" that we 
want to stay away from. 

• Back-propagate it for the appropriate time scale, M = [J <f>-t(N). The 

t£[0,T] 

counterpart to Fq should be chosen so that it does not enter M. M is not 
a small set; but we only need that cutting it off makes a small difference 
to Fq, not that itself it is small. Ffj(l — Xm( x , &)) is a starting point for 
F§; a simpler cutoff might be preferable. Certainly though, if ||Fq(1 — 
Xm(x, k))\\Li is not small, then this approach does not apply. 

• The previous step is only one constraint about at least how much we have 
to cut-off. If we cut-off any less (i.e. exclude a smaller set), we will enter 
too close to the singularity and the estimates will fail. An other constraint 
comes from having several H 2 bounds to check for our various approximate 
initial data. This means that if we exclude too small sets (e.g. too thin 
strips etc), the derivatives of the approximate data will be too large, and 
the estimates also fail. This is what makes us exclude a strip even when 
= in theorem |2.1[ where the previous step would be fulfilled by taking 
N to be a line. 

In other words, there are two different origins for constraints in this construction: 
on the one hand, we want to avoid the low-regularity region. The more we cut-off 
the better. On the other hand, we want the cut-offs to introduce small errors, 
including in various derivative-norms. This is a double sided constraint, as cutting 
off either a too thick or too thin strip will fail here. Of course the norms we can 
work in are constrained by the available conservation laws for the equations we 
work with, i.e. apparently have to be based on L 2 , and A' as long as we keep track 
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successfully of positivity. The fact that A' scales like L 1 in concentrating data is 
crucial in allowing wider strips to be cut-off. 

The fact that we will have to exclude domains measured in —(log e)~ l (to some 
power) means that in general one can only work with slowly concentrating initial 
data. Moreover, the repulsive character of the singularity is used implicitly here, 
as it helps the set M not be too large, but, when projected on the hyperplane 
transverse to propagation, basically looks like the set N. It is not necessary per se, 
but a way to meaningfully control how much larger is the set M from N is needed. 
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